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ABSTRACT 

This paper deals with the computational analysis of 
buoyancy effects in the near field of an isothermal 
helium jet injected into quiescent ambient air 
environment. The transport equations of helium mass 
fraction coupled with the conservation equations of 
mixture mass and momentum were solved using a 
staggered grid finite volume method. Laminar, 
axisymmetric, unsteady flow conditions were 

considered for the analysis. An orthogonal system with 
non-uniform grids was used to capture the instability 
phenomena. Computations were performed for Earth 
gravity and during transition from Earth to different 
gravitational levels. The flow physics was described by 
simultaneous visualizations of velocity and 
concentration fields at Earth and microgravity 
conditions. Computed results were validated by 
comparing with experimental data substantiating that 
buoyancy induced global flow oscillations present in 
Earth gravity are absent in microgravity. The 
dependence of oscillation frequency and amplitude on 
gravitational forcing was presented to further quantify 
the buoyancy effects. 

Keywords: Buoyancy, Flow computations. Flow 
Instability. 


INTRODUCTION 

Low density gas jet flows are encountered in many 
engineering and technical applications such as 
smokestack discharge, fuel leak, space propulsion, 
aerosol spray, jets entering mixing chamber, ventilation 
jets and many other environmental issues relating to 
atmosphere and hydrosphere. Jets are termed as 

buoyant for Richardson number, Ri= gd(p., -p)/pV" >1.0 
and inertial for Ri « 1.0. The present study examines 
the near field flow structure of a buoyant helium jet 
injected into quiescent ambient air environment. 
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Several studies 1 " 6 have focused on the far field behavior 
of jets. However, the near field flow stability has direct 
influence on the entrainment and flow pattern in the far 
field. Self-excited oscillations in the near injector 
region of low-density gas jets have been studied 
extensively. Subbarao and Cantwell observed highly 
periodic oscillations in the near field of helium jets 
injected into a co-flow air environment for Ri between 
0.5 and 6.0. The Strouhal number, St=fd/V correlated 
with Richardson number for Ri>1.0 indicating 
buoyancy dependent instability mode. Similar 
experiments were conducted by Cetegen 8 ' 10 for 
axisymmetric and planar plumes, whereby buoyancy 
induced toroidal vortical structures contaminating the 
primary jet flow were observed. Agrawal and co- 
workers 11 " 12 conducted experiments to characterize the 
flow instability in terms of concentration measurements 
across the whole field of a helium jet injected into 
quiescent air. Experiments by Yep et al. 11 demonstrated 
that flow oscillations were absent in the microgravity 
environment of the 2.2 s drop tower. 

Few computational studies have also been a part of 
low-density jet instability literature. Mell et al. 1 ' 
performed numerical simulations of helium jet in air 
and found that the computed oscillation frequency 
matched with experimentally obtained value by Hamins 
et al 14 . Recently, Soteriou 15 revealed that the instability 
characteristics in buoyant planar plumes compared well 
with experiments, providing insight into the role played 
by viscous and buoyancy forces in establishing 
oscillating or non-oscillating modes. 

The above studies provided extensive details relating to 
flow instability, attributed to buoyancy dominating the 
flow field. Although the instability in high Richardson 
number helium jets is attributed to buoyancy, direct 
physical evidence with respect to concurrent 
visualizations of velocity and concentration fields is 
lacking in the present literature. 
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If the instability is indeed caused by buoyancy, the near 
field flow structure of the jet will change significantly 
when buoyancy is removed. Thus, quantitative data on 
the spatial and temporal evolutions of the instability, 
length and time scales of the oscillating mode and its 
effects on the mean flow are needed in Earth gravity 
and microgravity conditions to delineate gravitational 
effects in buoyant jets. 

Although microgravity experiments 11 have quantified 
buoyancy effects and provided details of the 
concentration field, results are devoid of the flow 
structure details. Simultaneous visualization of velocity 
and concentration fields is desired to fully understand 
the oscillating instability and consequent flow structure 
of the jet. Owing to the experimental limitations in 
obtaining such measurements, especially in 
microgravity, the present study utilizes computational 
analysis to simulate helium jet in air. Features of the jet 
are described in both Earth gravity and microgravity. 
In the following section, we provide details relating to 
the governing equations, boundary and initial 
conditions, and the computational procedure. Results 
obtained at Earth gravity and microgravity conditions 
are discussed subsequently. 


Table 1. Coefficients and source terms in Equation 2. 


<!> 

T* 


W 

M 

dp . d ( dw^\ d ( dv) 2 d ( dv\ d ( dw) 

-r— + r (p^-p) g+ — rp — |+— 1 rp— rp— +— rp — 

dz dz[^ dzj dr\^ dz ) 3|_&\ dry dzl^ dzJJ 

V 

M 

dp 2 v df dvv) df dv] 8 f dv] df dwi] 

f 

Db 

0 


Thermo-Physical Properties 


The density was computed from the incompressible 
ideal gas law as follows: 
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The dynamic viscosity was calculated using the Wilke’s 
mixture averaged formula 16 as follows: 


NUMERICAL FORMULATION 


Governing Equations 

The governing equations for laminar, unsteady and 
axisymmetric flow are expressed in the cylindrical 
coordinate system as follows: 


at dr dz dr\ 
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where p represents density; and w and v are, 
respectively, the axial and radial components of the 
velocity. Equation 1 represents the continuity equation 
and Equation 2 is the general transport equation, which 
reduces to momentum, or species conservation equation 
depending on the variable^. The terms T^and sAn 
the transport equation represent the diffusion coefficient 
and source term of the variable <j> . The coefficients and 
source terms appearing in Equation 2 are listed in Table 
1 . In this study, the temperature and pressure variations 
are considered to be negligible and the jet flow is 
assumed to be a binary fluid system. 
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where subscripts 1 and 2 represent helium and air, 
respectively. The binary diffusion coefficient of helium 
into air, Db was specified as a constant. 

Operating and Boundary Conditions 

Figure 1 shows the computational domain extending 
15d in the axial direction and 6d in the radial direction 
where d=3 1.8mm is the tube inside diameter. The 
upstream boundary was placed at a distance of Id from 
the jet exit plane. Symmetric boundary condition was 
imposed along the jet centerline. At the tube inlet, the 
velocity profile was fully developed and the helium 
mole fraction was unity. The jet Reynolds number was 
specified as 300 corresponding to Richardson number 
of 1.52. No slip and zero mass diffusion conditions 
were imposed at the tube walls. At the exit boundary, 
pressure outlet condition was implemented instead of 


2 

American Institute of Aeronautics and Astronautics 




outflow conditions to account for the backflow. Figure 
1 shows pressure conditions imposed at other 
boundaries. The operating pressure and temperature in 
the domain were specified, respectively, as 1 atm and 
300 K. 



Computational Details 

The grid was generated by splitting the computational 
domain into five sub regions. An orthogonal non- 
uniform grid system was used. The grid elements in the 
radial direction were concentrated in the oscillatory 
zone. The grid size was increased gradually in the axial 
direction. The governing equations were solved 
sequentially using a segregated approach 17 . Since the 
governing equations are coupled, several iterations of 
the solution loop were performed before a converged 
solution was obtained. The pressure velocity coupling 
was achieved using the SIMPLEC algorithm. 


RESULTS AND DISCUSSION 

Numerical simulations were performed for the 
isothermal helium jet injected into the ambient air 
environment with the objective of isolating 
gravitational effects in the near field. First, grid 
sensitivity analysis was performed using different sizes 
of grids and the solution was found to be independent 
of the grid structure. Based on the analysis, a 260 x 95 
grid was chosen for computations. In the following 
sections, we provide details pertaining to the self 
excited instability observed in the helium jet at Earth 


gravity conditions. Subsequently, the near field flow 
physics is explained as the jet undergoes change to 
microgravity conditions. Further insight is gained from 
the parametric study by varying the gravitational level. 
The change in the buoyant structure of the jet is 
illustrated from vortex convection velocity, and flow 
oscillation frequency and amplitude plots. 

Self-Excited Oscillatory Structure in Earth Gravity 

Periodic oscillations occur in the near field of the 
helium jet at Earth gravity conditions owing to 
buoyancy induced acceleration of the jet fluid. To 
illustrate this phenomenon. Figure 2 shows a sequence 
(l/73s apart) of helium mole percentage contours 
superimposed with velocity vectors. The jet boundary 
is demarcated by the 60-70% concentration level and 
the vortex core is represented by a black dot. At t — 0.0s, 
buoyancy accelerates the jet core, which contracts to 
conserve the mass. The net result is entrainment of 
ambient fluid illustrated by the flow recirculation 
representing the vortex. In the next plot at t=l/73s, the 
vortex is characterized by larger recirculation region 
accompanied with greater contraction of the jet 
boundary, evident from the indentation in the 
concentration field. Subsequently, between t=2/73s and 
4/73s, the vortex convects downstream while growing 
in size to further contaminate the jet core. 

The vortex formation is accompanied by momentum 
transfer from the jet core to the vortex. This feature is 
evident in Figure 2 wherein vectors downstream of the 
vortex have higher magnitude than those upstream 
before interacting with the jet core. As the straining 
effect of the vortex diminishes in the upstream region, 
the jet expands at the exit inducting buoyant jet fluid to 
initiate another vortex at t=5/73s to repeat the 
oscillation cycle. The above results illustrate self 
excited periodicity in the jet at a frequency of 14.6 Hz, 
indicating both qualitative and quantitative agreement 
with experiments 11 . In the following sub-section, our 
study is focused on delineating buoyancy effects by 
simulating the jet flow in the absence of gravity. 

Transition from Earth to Microgravity Conditions 

Computations were performed to visualize the jet flow 
during change from Earth to microgravity conditions 
and to depict how the jet flow structure adjusts itself in 
the absence of buoyancy driven forces. As a first step, 
the simulation results were compared with 
experiments 11 to validate the computational model. 
Figure 3 shows a series of computational-experimental 
helium mole fraction contour plots during transition 
from Earth gravity to microgravity. Note that the 
experimental data were taken from Yep et al 11 . 
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Figure 2. Contours of helium mole percentage superimposed by velocity vectors during an oscillation cycle 
(left to right) in Earth gravity (a) t=0.0 s; (b) t=l/73 s; (c) t=2/73 s; (d) t=3/73 s; (e) t=4/73 s; (f) t=5/73 s. 


At t=0.0s, we observe a vortex located at z/d= 0.7, 
indicated by the indentation in the concentration level 
of 50%. With the initiation of microgravity, between 
t=l/60s and 2/60s, the vortex convects downstream as 
the jet expands in the near exit region, similar to 
experiments 11 . At t=3/60s, the jet has expanded at the 
exit to r/d= 0.7. The vortex located at z/d=1.5 has 
weakened as seen from the diminishing nature of the 
indentation and decrease in the curvature of 
concentration profile. Subsequently, the jet widens 
gradually throughout the near field with weakening and 
diminishing of oscillations as steady microgravity 
conditions are reached. The helium mole fraction 
contours in microgravity are straight lines reminiscent 
of a non-buoyant jet. The above mentioned features 
replicate similar phenomena observed in experiments 11 
wherein the self-excited oscillations gradually subsided 
after microgravity condition was initiated. Judging from 
the plots in Figure 3, it was inferred that the 
computations accurately reproduced the experiments 
proving the model’s reliability in predicting the flow 
physics during change from Earth gravity to 
microgravity. 

Figure 4 shows the contours of helium mole percentage 
overlapped by velocity vectors during transition from 
Earth to microgravity to highlight the interaction of the 
flow field with the concentration field. At t=0.0s, a 
vortex with the core located at z/d=0.7, characterizes 
the entrainment of ambient fluid into the helium jet. At 


this instant, the jet width near the exit extends to 
r/d=0.6. After the initiation of microgravity, at t=2/60s, 
the jet has expanded at the exit to r/d=0.7. The vortex 
convects downstream with its core located at about 
z/d=1.8. Upstream of the vortex, we observe an 
irregular distribution of the flow vectors together with 
distortions around the vortex core region. This behavior 
is attributed to the flow adjusting to the sudden decay in 
buoyancy induced acceleration along the jet core 
region. In the next plot at t=4/60s, the weakened vortex 
moves away from the field of view as another vortex is 
initiated at the jet exit. Observe that the vortex core is 
located close to the 90% concentration level with the jet 
expanding to r/d=0.75. Though buoyancy forces are 
absent, we observe vortical motion near the exit 
possibly because of the residual effects. 

At t=8/60s, we observe that the recirculation and 
entrainment have decreased with the jet continually 
expanding near the exit. Similar features are exhibited 
at t=10/60s with flow approaching steady conditions, 
whereby entrainment is dominated by diffusion. 
Finally, at t=2.0s, the jet has assumed a steady 
columnar shape. Higher velocities persist in the jet core 
and the velocity decreases radially as the jet mixes with 
the ambient air by molecular diffusion. The vorticity 
magnitude (not shown) decreases when buoyancy was 
removed, signifying reduction in mixing of helium and 
air by entrainment. Further information about the flow 
structure is provided by the spatio-temporal plots. 
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Figure 3. Measured (1) and Computed (2) Helium mole fraction contours during transition from Earth to microgravity at Re = 300, d = 31.8 mm, and 
Ri = 1.52. (a) t=0.0 s; (b) t=l/60 s; (c) t=2/60 s; (d) t=3/60 s; (e) t=13/60 s; (f) t=2.0 s. Measurements are taken from Yep et al u . 
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Figure 4. Contours of helium mole percentage superimposed by velocity vectors during transition from Earth 
gravity to microgravity (left to right) for (a) t=0.0 s; (b) t=2/60 s; (c) t=4/60 s; (d) t=8/60 s; (e) t= 10/60 s; (f) 
t=2.0 s. Note that gravity was removed at t=0.0sec. 


Temporal Evolution of Axial Velocity and 
Concentration Fields 

Figure 5 shows time trace plots of helium mole 
percentage and axial velocity at selected axial and 
radial locations during the change from Earth to micro 
gravity conditions. Near the jet exit at z/d=0.15, the 
concentration at the jet center is constant at 100% 
throughout. In contrast, the axial velocity at the jet 
center fluctuates between 2.3m/s and 2.6m/s in Earth 
gravity. With the initiation of microgravity conditions, 
the axial velocity initially decreases and then reaches a 
steady value of 2.2m/s within 0.1s of the change in 
gravity. In the entrainment region at r/d=0.45, the 
helium mole percentage in Earth gravity varied between 
45% and 94%. During transition to microgravity, the 
concentration peaks to 100% within 0.08s. Afterwards, 
the flow oscillates until t=0.2s with gradual decrease in 
the amplitude. Steady state was reached at t — 0.3s with 
helium concentration reaching a value of 95%. 
Correspondingly, the axial velocity fluctuations in 
Earth gravity diminish in microgravity, reaching a 
steady axial velocity of 0.4m/s. Near the jet boundary 
at r/d=0.6, low amplitude oscillations in the 
concentration level were observed in Earth gravity. 
Helium mole percentage reached a peak value of 70% 
within 0.1s after change of gravity. Subsequently, low 
Amplitude concentration oscillations were observed till 
t — 0 . 4 s before a steady value of 52% was reached at t = 


1.0s. Axial velocity fluctuations at this radial location in 
Earth gravity show troughs with negative velocity 
indicating flow recirculation. With the initiation of 
microgravity, the flow stabilized and reached a near 
zero axial velocity. 

At a downstream location of z/d=2.0, the helium 
concentration at r/d=0.0 remained constant throughout 
at 100% level. However, high amplitude oscillations 
are evident in axial velocity ranging between 2.4m/s 
and 4.6m/s in Earth gravity. With the initiation of 
microgravity, the fluctuations diminish within t — 0. 5s 
and reach steady conditions with an axial velocity of 
2.1m/s. At r/d=0.45, high amplitude concentration 
oscillations are observed in Earth gravity with helium 
concentration reaching peak levels of 95% during 
transition to microgravity and then reaching steady state 
within 0.6s. At r/d=0.6, helium reaches peak 
concentration level of 60% during transition to 
microgravity and reaches steady condition with 56% 
concentration after 0.6s. The velocity fluctuations at 
the above radial locations depict a similar flow behavior 
as already observed at the upstream location. 

The above observed features are in excellent agreement 
with time traces of concentration fields obtained by Yep 
et al. 11 using the 2.2 sec drop tower. The jet expanded 
in microgravity and flow oscillations in concentration 
field existed for a brief period before attaining steady 
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Figure 5. Time traces of helium mole percentage and axial velocity during transition from Earth gravity to 
microgravity at a) z/d=0.15, and b) z/d=2.0. 


state conditions. Observe that an upstream location 
reaches steady state prior to the downstream location. 
Coupling between the concentration and velocity 
fields is explained next using radial profiles near the 
jet exit (z/d=0.15) at Earth and microgravity 
conditions to focus on the buoyancy effects in the jet 
shear layer. 

Figures 6 and 7 depict the radial profdes of 
normalized axial velocity and helium mole 
percentage at z/d=0.15 for Earth and microgravity 
conditions. Note that the axial velocity has been 
normalized with the jet exit velocity at respective 
radial locations. In Earth gravity, the normalized 
axial velocity profile shows a gradual increase from 
the jet core towards the tube wall indicating greater 
acceleration in the outer region of the jet. The peak 
normalized axial velocity corresponds to a value of 
1.2 at r/d=0.4, a location coinciding with the sudden 
drop in concentration. Thus, we observe that 
buoyancy causes large axial velocity gradients and 
unstable density stratification near the jet shear layer, 
leading to the entrainment of ambient fluid into the 
primary jet to produce the oscillatory flow. 


In microgravity, a flat normalized axial velocity 
profile with a minor dip near the tube wall is noticed 
(Ref. figure 7). The concentration changes gradually 
as helium diffuses into the ambient environment. 
Results show that the lack of buoyant acceleration in 
the jet shear produces a non-oscillatory flow. The 
next section provides quantitative details relating to 
the vortex movement and frequency and amplitude of 
flow oscillations with changes in the gravitational 
forcing. 

Flow Structure: Dependence on Gravitational Forcing 

In this section, the analysis has been limited only to 
the vortex movement and the frequency and 
amplitude of flow oscillations as functions of the 
gravitational forcing. Figure 8 illustrates the 
convection velocity of the vortex core at different 
gravitational levels. For Martian gravity ( g/go — 0 .37), 
we observe that the vortex convected at a uniform 
acceleration with peak velocity of 0.9 m/s at z/d=3.0. 
Further, in Earth gravity (g/go=1.0), the vortex 
convection velocity increased at all the axial stations 
with a peak velocity of 1.0 m/s. 
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Figure 6. Helium mole percentage and normalized Figure 8. Vortex convection velocity at various 

axial velocity profiles in Earth gravity at z/d=0.15. gravitational levels. 



r/d Normalized gravitational acceleration, (g/go) 


Figure 7. Helium mole percentage and normalized Figure 9. Oscillation frequency Vs normalized 

axial velocity profiles in microgravity at z/d=0.15. gravitational acceleration. 


In Jupiter gravity (g/g 0 =2.33), the vortex convects 
with higher acceleration throughout the near field as 
shown by the steep rise of the curve, with peak 
velocity reaching about 2.0m/s. Notice that with 
increase in buoyancy, the vortex convection velocity 
increases throughout the near field. Accordingly, the 
entrainment rate, which is a function of the vortex 
convection velocity, would also increase at higher 
gravitational levels. Higher entrainment rate leads to 
enhanced mixing and contamination of the jet core in 
the downstream region. 

Figures 9 and 10 show plots explaining the variation 


of frequency and amplitude of flow oscillation with 
the gravitational level. The best fit resulted in the 
relationship between oscillation frequency and 
gravitational acceleration as f = 9.38(g/go) 02 . The 
normalized amplitude of the jet oscillatory structure 
in Figure 10 shows dependence on the gravitational 
level as well. The amplitude refers to that of the 
concentration wave and normalization is done using 
amplitude at Earth gravity. An increase in amplitude 
with the rise in gravitational level suggests that 
buoyancy enhances wrapping of the ambient fluid 
into the jet flow field to affect the entrainment and 
mixing characteristics in the oscillating shear layer. 


8 


American Institute of Aeronautics and Astronautics 





2.5 


NOMENCLATURE 



Normalized gravitational acceleration, (g/go) 


Figure 10. Normalized concentration amplitude 
Vs normalized gravitational acceleration at 
z/d=1.0 and r/d=0.45. 


CONCLUSIONS 


Numerical simulations of a helium jet injected into 
quiescent air environment were performed to 
describe the effects of buoyancy in the jet flow field. 
The interaction of the velocity and the concentration 
fields was presented by simultaneous visualizations 
at Earth and microgravity conditions. Results 
indicated gravity induced global oscillations 
occurring in the near field of the jet at Earth gravity. 
Further, the gradual weakening and diminishing of 
oscillations during transition from Earth to 
microgravity confirmed that flow oscillations are 
buoyancy induced. Temporal plots of velocity and 
concentration fields depict that oscillations existed 
for a brief period during transition from Earth to 
microgravity. The jet width increased in microgravity 
and the upstream locations of the jet reached steady 
state prior to the downstream locations. Radial 
profiles of normalized axial velocity and 
concentration fields near the jet exit illustrated the jet 
diffusion and entrainment processes characterizing 
the instability. It was found that the flow oscillation 
frequency and amplitude increase with the rise in the 
gravitational level. 
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g Accelaration due to gravity (m/s 2 ) 

M w Molecular Weight (Kg/Kmole) 
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Subscripts 
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he Helium 

i Species i 

j Species j 

co Ambient fluid 
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